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Description 



[Interpretation Methods for NMR 
Diffusion-T2 maps] 

Cross Reference to Related Applications 

[0001] This invention claims priority pursuant to 35 U.S.C. §119 
of U.S. Provisional Patent Application Serial No. 
60/450,412, filed on February 27, 2003. This Provisional 
Application is hereby incorporated by reference in its en- 
tirety. 

Background of Invention 

[0002] Both water and hydrocarbons in earth formations produce 
detectable NMR signals. It is desirable that the signals 
from water and hydrocarbons be separable so that hydro- 
carbon-bearing zones may be identified. However, it is 
not always easy to distinguish which signals are from wa- 
ter and which are from hydrocarbons. Various methods 
have been proposed to separately identify water and hy- 
drocarbon signals. 

[0003] The differential spectrum (DSM) and shifted spectrum 



(SSM) methods proposed by Akkurt et al in "NMR Logging 
of Natural Gas Reservoirs " Paper N. Transactions of the 
Society of Professional Well Log Analysts (SPWLA) Annual 
Logging Symposium, 1995, compare 1^ distributions de- 
rived from two Carr-Purcell-Meiboom-Gill (CPMG) mea- 
surements performed with different polarization times 
(DSM) or echo-spacings (SSM). A modification to these 
methods, known as time domain analysis (TDA), was later 
introduced by Prammer et al. in "Lithology-lndependent 
Gas Detection by Gradient-NMR Logging, "SPE paper 
30562, 1995. In TDA, "difference"data are computed di- 
rectly in the time domain by subtracting one set of the 
measured amplitudes from the other. The difference 
dataset is then assumed to contain only light oil and/or 
gas. In TDA, relative contributions from light oil or gas are 
derived by performing a linear least squares analysis of 
the difference data using assumed NMR responses for 
these fluids. Both DSM and TDA assume that the water 
signal has substantially shorter^ relaxation times than 
those of the hydrocarbons. This assumption is not always 
valid, however. Most notably, this assumption fails in for- 
mations where there are large pores or where the hydro- 
carbon is of intermediate or high viscosity. The SSM 



method and its successor, the enhanced diffusion method 
(EDM) proposed by Akkurtef. al in "Enhanced Diffusion : 
Expanding the Range of NMR Direct Hydrocarbon Typing 
Applications", Paper GG. Transactions of the Society of 
Professional Well Log Analysts (SPWLA) Annual Logging 
Symposium, 1998, separate gas, oil and water contribu- 
tions based on changes in the distributions that result 
from changes in the echo spacing of CPMG measure- 
ments. The methods are applicable in a limited range of 
circumstances and the accuracy of the result is signifi- 
cantly compromised by incomplete separation of water 
and hydrocarbon signals in the domain. Moreover, 
these methods are designed to function with CPMG se- 
quences. However, with the diffusion-based methods, 
CPMG pulse sequences provide poor signal to noise ratios 
due to the reduced number of echoes that can be mea- 
sured. A strategy for combining and selecting these dif- 
ferent NMR methods has been described recently by 
Coates et al. in U.S. Patent No. 6,366,087 Bl. 
[0004] The diffusion-editing (DE) pulse sequence by Hurlimann et 
al. provides a different approach. See M.D. Hurlimann et 
al . , "Diffusion-Editing: New NMR Measurement of Saturation and 
Pore Geometry "paper presented at the 2002 Annual Meet- 



ing of the Society of Professional Well Log Analysts, Osio, 
Japan, June 2-5; see also, U.S. Patent No. 6,570,382, filed 
on November 28, 2000, by Hurlimann. This patent is as- 
signed to the same assignee as the present invention and 
is hereby incorporated by reference. DE pulse sequences 
are similar to the CPMG sequences except that the initial 
two echoes are acquired with longer echo spacings and 
the third and subsequent echoes are acquired with shorter 
echo spacings. In DE pulse sequences, diffusion informa- 
tion is encoded during the acquisition of the first two 
echoes, whereas the third and subsequent echoes provide 
bulk and surface relaxation time information with rela- 
tively little attenuation of the signal by diffusion. Using a 
conventional CPMG sequence to encode the diffusion in- 
formation requires a long inter-echo spacing, which re- 
sults in poor bulk and surface relaxation time information 
because diffusion decay attenuates the signal after rela- 
tively few echoes. Consequently, a suite of data acquired 
with DE sequences provides better diffusion information 
and signal-to-noise ratio in the spin-echo data, as com- 
pared to an analogous suite acquired with CPMG se- 
quences. Therefore, DE sequences can provide more ac- 
curate and robust computations of brine and oil T distri- 



butions than CPMG sequences. 

[0005] | n addition to DE sequences, specialized interpretation 
methods have been developed for NMR data in order to 
further enhance hydrocarbon detection. These methods 
typically apply forward modeling to suites of NMR data 
acquired with different parameters. The suite of NMR data 
are typically acquired with different echo spacings (TE) or 
polarization times (WT), and sometimes acquired with dif- 
ferent magnetic field gradients (G). DE sequences are one 
example of such data acquisition. Two exemplary meth- 
ods include: the MACNMR proposed by Slijkerman etal, 
SPE paper 56768, "Processing of Multi-Acquisition NMR 
Data", 1999, and the Magnetic Resonance Fluid character- 
ization (MRF) methoddisclosed in U.S. Patent No. 
6,229,308 Bl issued to Freedman and assigned to the as- 
signee of the present invention ("the Freedman patent"). 
The Freedman patent is hereby incorporated by reference. 

[0006] The MRF method is capable of obtaining separate oil and 
water T distributions. This method uses a Constituent 

2 

Viscosity Model (CVM), which relates relaxation time and 
diffusion rates to constituent viscosities whose geometric 
mean is identical to the macroscopic fluid viscosity. With 
the MRF method, estimates for water and hydrocarbon 



volumes are obtained by applying a forward model to 
simulate the NMR responses to a suite of NMR measure- 
ments acquired with different parameters. Specifically, 
The MRF technique is based on established physical laws 
which are calibrated empirically to account for the down- 
hole fluid NMR responses. By using realistic fluid models, 
MRF aims to minimize the number of adjustable parame- 
ters to be compatible with the information content of typ- 
ical NMR log data. Since the model parameters are by de- 
sign related to the individual fluid volumes and proper- 
ties, determination of the parameter values (i.e. data- 
fitting) leads directly to estimates for petrophysical quan- 
tities of interest. 
[0007] The forward-model approach relies on the validity of the 
fluid models employed. In "non-ideal" situations where 
fluid NMR responses deviate from the model behavior 
(oil-wet rocks, restricted diffusion), these techniques may 
lead to erroneous answers. In some circumstances, "non- 
ideal'Vesponses may be identified by poor fit-quality, in 
which case the fluid models can be adjusted by modifying 
the appropriate model parameter. However, it may not be 
obvious which element of the fluid model should be mod- 
ified and what modification is needed. 



[0008] Another approach developed by Schlumberger, based on a 
maximum entropy principle (MEP), consists of a general 
model-independent method to analyze complex fluids 
data acquired with NMR logging instruments and present 
the results in a visually attractive and easy-to-understand 
format, hereby referred to as Diffusion-Relaxation maps, 
or D-T2 maps. These maps have been used to understand 
cases where model-based analysis gives unsatisfactory 
results because of deviations of NMR properties from the 
"ideal" behavior assumed in the models. These situations 
can arise due to anomalous fluid/rock interactions such as 
restricted diffusion, mixed-wettability and internal gradi- 
ents. Deviations from the default properties have also 
been observed for certain crude oils, leading to inaccurate 
predictions in the model analysis. Through the use of D- 
T2 maps, the MEP approach provides a simple graphical 
representation of the data that can be used to identify 
fluid responses in all environments. Diffusion-Relaxation 
maps are further described in commonly assigned U.S. 
Patent Nos. 6,570,382 and 6,462,542. 

[0009] while these prior art methods are useful in predicting the 
presence of hydrocarbons in the formations, it is desirable 
to have simpler methods that can predict the presence of 



hydrocarbons in the formations from NMR data and are 
generally applicable to NMR data acquired with different 
pulse sequences. Furthermore, while two and three di- 
mensional visualization has been developed to obtain pri- 
marily qualitative information, it is desirable to have 
quantitative interpretation techniques that can provide ac- 
curate fluid-characterization results. 
Summary of Invention 

[0010] According to one aspect of the disclosed subject matter a 
method is described for interpretation of multi-di- 
mensional nuclear magnetic resonance data taken on a 
sample of an earth formation. Specifically, a set of NMR 
data is acquired for a fluid sample located either in a 
borehole or in a laboratory environment. From the set of 
NMR data, a multi-dimensional distribution is calculated 
using a mathematical inversion that is independent of 
prior knowledge of fluid sample properties. The multi- 
dimensional distribution is graphically displayed on a 
multi-dimensional map. Each fluid instance or artifact vis- 
ible on the graph is identified as representing a probable 
existence of a detected fluid. One or more quantitative 
formation evaluation answers for one or more fluid in- 
stances are computed based on the multi-dimensional 



distribution associated with the respective fluid instance. 

[0011] According to another aspect, quantitative formation eval- 
uation answers are determined from the multi-di- 
mensional distribution of NMR data by initially determin- 
ing a set of model parameters which represent aspects of 
the multi-dimensional distribution. A model dependent 
inversion is then applied to compute the fluid properties. 

[0012] According to another aspect, quantitative formation eval- 
uation answers are determined from the multi-di- 
mensional distribution of NMR data through a point- 
and-click approach. One or more fluid artifacts are se- 
lected from a multi-dimensional map of the NMR data us- 
ing a computer mouse or an automatic edge selection ap- 
plication. The amplitude is integrated over the selected 
region to determine properties of the fluid associated with 
the selected region. 

[0013] According to another aspect, quantitative formation eval- 
uation answers are determined from the multi-di- 
mensional distribution of NMR data by determining a 
mean diffusion value across a region of a diffusion-T2 re- 
laxation distribution. The mean diffusion is used to deter- 
mine properties of the fluid associated with the selected 
region. 



Brief Description of Drawings 

[0014] FIG. 1 is a diagram of an exemplary downhole nuclear 

magnetic resonance data acquisition system. 
[0015] FIG. 2 is a more detailed diagram of the system of FIG. 1. 

[0016] FIG. 3 is a multi-dimensional map or graph for displaying 
NMR data. 

[0017] FIG. 4 is a NMR data map showing non-ideal effects on 
the NMR data. 

[0018] FIG. 5 is another NMR data map showing non-ideal effects 
on the NMR data. 

[0019] FIG. 6 is another NMR data map showing non-ideal effects 
on the NMR data. 

[0020] FIG. 7 is another NMR data map showing non-ideal effects 
on the NMR data. 

[0021] FIG. 8 is two dimensional graph according to a prior art 
inversion method. 

[0022] FIG. 9. is a NMR data map set showing correction of non- 
ideal effects according to one aspect of the disclosed sub- 
ject matter. 

[0023] FIG. 10 illustrate a correction of a prior art inversion 
method according to the NMR data map set of FIG. 9. 
[0024] FIG. 11 is two dimensional graph according to a prior art 



inversion method. 
[0025] FIG. 12. is a NMR data map set showing correction of 

non-ideal effects according to one aspect of the disclosed 

subject matter. 
[0026] FIG. 13 illustrate a correction of a prior art inversion 

method according to the NMR data map set of FIG. 12. 
[0027] FIG. 14 is a set of NMR maps and graphs according to one 

formation evaluation method of the disclosed subject 

matter. 

[0028] FIG. 15 is a flow diagram of the method illustrated in FIG. 
14. 

[0029] FIG. 16 is a multi-dimensional map or graph of another 
formation evaluation method of the disclosed subject 
matter. 

[0030] FIG. 17 is a set of NMR maps and graphs according to the 

method illustrated in FIG. 16. 
[0031] FIG. 18 is a flow diagram of the method illustrated in FIGS. 

16 and 17. 

[0032] FIG. 19 is another set of NMR maps and graphs of a for- 
mation evaluation method of the disclosed subject matter. 

[0033] FIG. 20 is a set of corrected NMR graphs according to the 
method illustrated in FIG. 19. 

[0034] FIG. 21 is a screen shot illustrating exemplary formation 



evaluation answers provided by the disclosed subject 

matter. 
Detailed Description 



[0035] The disclosed subject matter describes quantitative meth- 
ods to interpret two-dimensional nuclear magnetic reso- 
nance (NMR) maps derived from common NMR formation 
evaluation measurements. Although other values may be 
used, a preferred embodiment is primarily discussed 
herein based on diffusion vs. (D-T2) maps. According 
to the present invention, D-T2 maps can be used to help 
in the selection of parameters for application to existing 
model-based inversion codes. Further, complete petro- 
physical answers (porosity, permeability, fluids volumes, 
saturations, oil viscosity etc.) can be derived directly from 
the D-T2 maps. To take advantage of the visual appeal of 
the maps, the proposed methods are interactive and, ac- 
cording to an embodiment, consist of sequential point- 
and-click procedures. 

[0036] Acquisition of NMR measurements according to embodi- 
ments of the invention may be accomplished with various 
methods of NMR measurements known in the art. For ex- 
ample, the measurements may be performed in a labora- 
tory using a sample removed from an earth formation. Al- 



ternatively, the NMR measurements may be performed in 
a logging operation using a wireline tool, a logging- 
while-drilling or measurement-while-drilling tool, or a 
formation tester. FIG. 1 illustrates a schematic of an NMR 
logging system. In FIG. 1, a NMR logging tool 30 for in- 
vestigating earth formations 31 traversed by a borehole 
32 is shown. The NMR logging device 30 is suspended in 
the borehole 32 on an armored cable 33, the length of 
which substantially determines the relative axial depth of 
the device 30. The cable length is controlled by suitable 
means at the surface such as a drum and winch mecha- 
nism 8. Surface equipment 7 can be of conventional type 
and can include a processor subsystem which communi- 
cates with downhole equipment including NMR logging 
device 30. 

[0037] The NMR logging device 30 can be any suitable nuclear 
magnetic resonance logging device; it may be one for use 
in wireline logging applications as shown in FIG. 3, or one 
that can be used in logging-while-drilling (LWD) or mea- 
surement-while-drilling (MWD) applications. In addition, 
the NMR logging device 30 may be part of any formation 
tester known in the art, such as that sold under the trade 
name of MDT™ by Schlumberger Technology Corporation 



(Houston, TX). The NMR logging device 30 typically in- 
cludes a means for producing a static magnetic field in 
the formations, and a radio frequency (RF) antenna means 
for producing pulses of magnetic field in the formations 
and for receiving the spin echoes from the formations. 
The means for producing a static magnetic field may 
comprise a permanent magnet or magnet array, and the 
RF antenna means for producing pulses of magnetic field 
and receiving spin echoes from the formations may com- 
prise one or more RF antennas. 
[0038] FIG. 2 illustrates a schematic of some of the components 
of one type of NMR logging device 30. FIG. 4 shows a first 
centralized magnet or magnet array 36 and an RF antenna 
37, which may be a suitably oriented coil or coils. FIG. 2 
also illustrates a general representation of closely-spaced 
cylindrical thin shells, 38-1, 38-2. . .38-N, that can be 
frequency selected in a multi-frequency logging opera- 
tion. One such device is disclosed in U.S. Patent No. 
4,710,713. In FIG. 2, another magnet or magnet array 39 
is shown. Magnet array 39 may be used to pre-polarize 
the earth formation ahead of the investigation region as 
the logging device 30 is raised in the borehole in the di- 
rection of arrow Z. Examples of such devices are disclosed 



in U.S. Patent Nos. 5,055,788 and 3,597,681. 

[0039] Turning now to FIG. 3, shown is an exemplary D-T2 map 
with NMR spin echo data presented as amplitudes versus 
diffusion (D) and relaxation (T2). The map shown in the 
left panel is a three-axis perspective view. The right panel 
provides a more practical representation of D-T2 map as 
a two-axis map. However, it should be noted that the dis- 
closed methods may be applied a dataset having any 
number of dimensions, 2-D, 3-D, 4-D, etc. Furthermore, 
it should be noted that although D-T2 maps are discussed 
herein for exemplary purposes, the disclosed methods can 
be equally as effective in obtaining quantitative formation 
evaluation answers based on many other combinations of 
NMR data properties (D, Tl, T2, T1/T2, etc.). 

[0040] | n the context of the two-axis D-T2 map, the diffusion 
amplitude is represented according to a color-coding 
scheme. The differences of diffusion properties among 
gas, water and various viscosity oils are captured by the 
D-T2 map and shown as separate and distinct peaks. 
Specifically, the color grouping at A, also herein referred 
to as an artifact or a fluid instance, represents the proba- 
ble detection of a first fluid. Similarly, the lighter color 
groupings or fluid instances at B, C and D also represent 



the probable detection of three additional fluids. The the- 
oretical responses of water, oil, and gas are overlaid on 
the maps to help the interpretation. Thus, for grouping or 
instance A, it is likely that the fluid is gas because its peak 
lies near the theoretical gas diffusion value. For groupings 
B and C, it is likely the fluids are varying viscosities or 
phases of oil, lying along the theoretical oil diffusion line. 
Finally, it is likely that grouping or instance D is water 
subject to restricted diffusion (discussed below). 

[0041] According to one embodiment, the cross-plot of FIG. 3 is 
model-independent. This means that no predetermined 
diffusion values or limits are imposed in the calculation to 
obtain NMR measurements from the spin echo data. Al- 
though others exist, one example of a model-inde- 
pendent calculation is the previously mentioned MEP ap- 
proach. Returning to FIG. 3, in order to assist interpreta- 
tion of the NMR data, an overlay of theoretical responses 
of water, oil and gas is helpful. 

[0042] As mentioned, previous attempts to determine quantita- 
tive formation evaluation answers have been based on the 
model predetermined values of diffusion and T2 relax- 
ation parameters of the fluids. In addition, it has been re- 
quired to select a fluid model based upon a best guess as 



to which fluids will be detected. Needless to say, any inac- 
curacies in the initial estimates of the fluid model and the 
fluids parameters create inaccuracies in the final answers. 
FIG. 4 graphically illustrates some exceptions to common 
ideal models: 1) Internal gradient pulls data in the up- 
arrow direction (such that water could be mistaken for 
gas), 2) Restricted diffusion pulls data in the other direc- 
tion as indicated by the down arrow (such that water could 
be mistaken for oil), 3) Mixed wettability pulls the oil data 
to the left (such that the oil volume/saturation are com- 
puted too low), and 4) high GOR (Gas Oil Ratio)cause oils 
to shift along the North-East trend (such that oils could be 
mistaken for water or gas). 

[0043] For example, FIG. 5 illustrates an effect of internal gradi- 
ent effect on a NMR fluid response. The internal gradient 
adds to the tool gradient and hence, diffusion is higher 
than expected. The sample is a clay-rich sandstone in a 
known water-bearing zone of a well drilled with water- 
base mud. As such, the calculated diffusion values appear 
higher than the expected result, indicated by the water 
diffusion line overlay. 

[0044] Another example is shown in FIG. 6, which illustrates an 
effect of restricted diffusion in a known carbonate domi- 



nated formation. The free water undergoes unrestricted 
diffusion in large pores and agrees with the theoretical re- 
sponse. The bound water is trapped in smaller pores and 
therefore experiences restricted diffusion. The result is 
one peak, representing the free water, at the expected 
water diffusion overlay and a second peak, representing 
the bound water, stretched across a range of diffusion 
values, even crossing into values indicative of hydrocar- 
bon presence. 

[0045] yet another example, shown in FIG. 7, illustrates gas in 

high-permeability sandstone reservoirs in the Middle East. 
The large pore sizes cause the water to relax close to its 
bulk value overlapping the gas signal, and resulting in a 
smeared peak at a long T2. 

[0046] According to an embodiment of the present invention, 
these D-T2 maps are generated by model-independent 
inversion codes such as MEP. These model independent 
inversions do not require a priori input of the fluids diffu- 
sion or T2 values. According to one embodiment, the 
model independent inversions do not require any a priori 
knowledge of fluid properties nor what fluids are present. 
From these inversions, D-T2 maps or graphs are gener- 
ated to display the resultant multi-dimensional NMR data 



across two, three or more axes in an easily readable form. 
[0047] According to one embodiment, the D-T2 maps are used 
to improve the results of model-dependent inversion such 
as the MRF analysis. Specifically, MEP, and other model in- 
dependent derived D-T2 maps provide an unbiased rep- 
resentation of the NMR data. As such, an overlay of the 
theoretical responses of the three most encountered flu- 
ids, i.e., water, oil, and gas indicates whether the re- 
sponses need adjustments for use in model-based inver- 
sion. 

[0048] For example, the MRF model (for oil, gas and water) states 
that the water and gas diffusion constants are indepen- 
dent of T2, and depend on temperature rand pressure P 
(for gas): 

DJT2) = DJT) (1) 

D t (72) =D i (T,P) (2) 

For oil, the diffusion constant is linearly proportional to 

T2, 

DJTT) = X x T2 (3) 

[0049] it follows from Equations (1) - (3) that two horizontal lines 
(i.e. at constant D values) representing the theoretical re- 
sponses of water and gas, and a diagonal line represent- 
ing the theoretical response of oil can be overlaid on a D- 



T2 map. Deviations from the ideal fluid responses will be 
evident in the maps as signals located away from the 
overlay lines. Once known, these deviations are applied to 
the model-based inversion. The result from the inversion 
provides an answer with improved accuracy based on the 
observations from the model-independent D-T2 maps. 
[0050] FIGS. 8-10 illustrate the model-dependent inversion cor- 
rection in a water based drilling fluid, cesium formate, and 
FIGS. 11-13 illustrate the correction in a primarily crude 
oil sample. FIG. 8 shows an exemplary model dependent 
based analysis, such as the MRF approach, for a cesium 
formate solution. FIG. 8 shows the model-based analysis 
having an oil peak 200 and a water peak 201. The model 
based analysis using the default water diffusion constant 
indicates mostly oil for cesium formate drilling fluid. FIG 9 
shows the D-T2 map on the left for the same data having 
the fluid instance or artifact 206 lying along the default oil 
diffusion line 204. Specifically, the D-T2 map on the left 
shows that the detected Cs formate fluid, the instance 
206, has a much lower diffusion constant than pure water 
(the default), indicated by horizontal line 202. Thus it can 
be determined that the ideal water diffusion 202 is not 
correct for the Cs formate sample. Thus, in order to com- 



pute correct fluid saturations the effective "water" (Cs for- 
mate to be precise) diffusion constant must be reduced to 
40% of its default value, represented by horizontal line 
208. This corrected diffusion value is then programmed 
into the model based inversion. FIG. 10 shows the repro- 
cessed MRF results obtained using a water diffusion con- 
stant reduced to 40% of its original value. The MRF results 
now correctly indicate predominantly water, shown at 
peak 210, having a very low oil signal, shown at peak 212. 
As such, the associated formation evaluation determina- 
tions such as water volume and saturation can be more 
accurately computed. 
[0051] FIG. 11 shows a model dependent based analysis of a 

known crude oil sample. The results indicate the presence 
of oil, shown by peak 220, but also incorrectly indicates 
significant water saturation, indicated by peak 221. FIG. 
12 shows the D-T2 map for the same data. On the left, 
the D-T2 map is shown with the default fluid overlay 
lines, the water diffusion represented by horizontal line 
222 and the oil diffusion represented by line 224. The 
signal or instance 226 appears above the oil D-T2 corre- 
lation line 224 indicating that this oil has an unusually 
high diffusion to T2 ratio. Thus, it can be determined that 



the existing oil D-T2 line 224 is incorrect and requires 
adjustment. The same map is plotted on the right with a 
new oil D-T2 correlation line 228 that bisects the main 
signal peak. The MRF analysis is then run using an in- 
creased 

value (see Eq (3)). The results of the reprocessed data are 
shown in FIG 13. As expected, the MRF analysis now cor- 
rectly predicts predominantly oil at peak 232, while indi- 
cating a very low water presence at curve 230. 
[0052] | n addition to assisting model-dependent interpretation 
techniques as shown above, complete quantitative petro- 
physical answers can be derived directly from the two- 
dimensional maps. Specifically, according to one applica- 
tion of the invention, quantitative measurements of 
porosity, permeability, fluid volumes, saturations, oil vis- 
cosity and other quantities can be derived from D-T2 
maps. Compared to prior methods, additional interpreta- 
tion is needed to derive more than previously utilized 
qualitative information. According to the disclosed subject 
matter, two approaches, a point-and-click approach and a 
diffusion log mean approach, are used to obtain quantita- 
tive answers. 



[0053] According to one embodiment, a visual point-and-click 

approach is provided which allows a user to interact with a 
D-T2 map by focusing in on particular artifacts graphi- 
cally shown on a map. The signal amplitude, A, from a 
suite of NMR pulse sequences can be expressed as 

(4) 

where WT, TE, fare the wait time, echo spacing and time of 
the NMR pulse sequences, (i, jk) are the indices of the T2, 
D and T1/T2 distributions, f(ij,k) is the amplitude of the 
three-dimensional component in the T2, D, T1/T2 space, 
H(WT f TE, t i, J, k) is the kernel of that component, and □ is 
a noise term. 

[0054] prom the above equation (4), it can be seen that the D-T2 
map is a representation of the signal amplitudes inte- 
grated across the k dimension (T1/T2) of the MEP inver- 
sion result. Therefore, in the favorable case where the flu- 
ids distributions in the D-T2 space are well separated as 
seen in FIG. 14, their respective volumes can be directly 
obtained by integrating the amplitudes in the D and T2 
windows defined by each type of fluid. The respective sat- 
urations can then be derived by dividing the respective 
fluid volumes with the total fluid volume. The individual 



volumes must be corrected for hydrogen index, according 
to known methods, to give correct answers. 

[0055] FIG. 14 illustrates an example according to one embodi- 
ment of the point-n-click approach having four fluids ar- 
tifacts, 240, 242, 244 and 246, that are well separated in 
D-T2 space (top left pane). Note that the four fluids are 
not resolved in either T2 space (bottom left pane) or D 
space (top right pane). Specifically, the T2 and D space 
graphs indicate three instead of four fluids having separa- 
ble properties. In the bottom left 12 graph, peak 248 cor- 
responds to artifact 246 and peak 250 corresponds to ar- 
tifact or fluid instance 244. However, peak 252 cannot re- 
solve artifacts 240 and 242. In the top right diffusion dis- 
tribution, peak 254 corresponds to artifact 240 and peak 
256 corresponds to fluid instance 242. However, peak 
258 cannot resolve artifacts 244 and 246. 

[0056] According to prior methods, quantitative fluid answers 

could be obtained using T2 graphs. However, as shown, in 
some cases T2 maps cannot fully resolve multiple fluids 
having different diffusion properties. Known prior meth- 
ods have gone a step further to use the D-T2 map to 
evaluate the accuracy (qualitatively) of the answers derived 
from the T2 graphs. According to an embodiment of the 



present invention, a visual point-n-click method uses the 
D-T2 map to determine quantitatively the respective vol- 
umes of each of the four fluids by integrating the signal 
amplitudes along 12 and D dimensions in the windows 
defined by the rectangles 260, 262, 264 and 266. Note 
that the use of other shapes to delineate the map region 
of interest, such as polygons or circles, may be employed 
in a similar manner. In this way, the disclosed methods 
advance the state of the art in part by resolving multiple 
fluids having similar T2 distributions and determining 
quantitatively certain formation evaluation answers once 
the fluid artifacts are separately identified. 
[0057] According to one application of the point-and-click 

method, the interpreter can easily select the integration 
region of the D-T2 map using, for example, the computer 
mouse or a digitized pen. Selection of a region may also 
be performed automatically by a software algorithm, for 
example, based on a predetermined amplitude threshold. 
Note also that the interpretation of the fluid type is guided 
by overlaying the theoretical responses of gas, oil and wa- 
ter in the D-T2 map as seen in the top left pane. This step 
also may be performed automatically by a software appli- 
cation, for example, based on a proximity of a point of 



maximum amplitude to the theoretical gas, oil and water 
responses. 

[0058] Turning to FIG. 15, illustrated is a flow diagram of the 

point-n-click approach. Beginning at step 500, the point- 
and-click routine begins with a D-T2 map. As mentioned, 
this D-T2 map is preferably generated using a model- 
independent inversion approach. Continuing at step 502, 
the user is asked to enter the type of fluid model to be 
used in the calculation. Typically, this will be done by in- 
spection of the map. First, the user decides on the num- 
ber of fluids. For example, if two amplitude clusters 
(bright spots) are present, then a two-fluid approach is 
chosen. Second, each fluid type is interpreted with the 
help of the overlay of theoretical responses of water, oil 
and gas. At step 504, the user selects the fluid region us- 
ing a computer mouse. According to one embodiment, a 
box is displayed and manipulated to encompass substan- 
tially all the artifact on the D-T2 map. According to a em- 
bodiment, the final unselected region having a positive 
non-zero amplitude is accumulated and displayed to pro- 
vide an indication to the user as to whether any substan- 
tial portions had been missed. According to another em- 
bodiment, a computer or software application selects an 



extent at least partially surrounding of each fluid instance 
based on a threshold amplitude. 
[0059] Once the artifact has been selected, the fluid volume and 
saturation is calculated at step 506 by integrating over the 
selected region of the D-T2 map. Since the integration 
over the total map area gives total porosity, individual 
fluid saturations can be computed by dividing the fluid 
volumes with total porosity. It is appropriate to reiterate 
here that although the disclosed point-and-click method 
is discussed for exemplary purposes in the context of D- 
T2 maps, nearly any multi-dimensional map can be em- 
ployed to determine a quantity of interest. Continuing at 
step 510, once fluid volume and saturation quantities, or 
other base values, have been determined, auxiliary com- 
putations at 510 can be requested by the user or auto- 
matically by the algorithm. For example, at step 512, vis- 
cosity can be determined by computing the mean 12 in 
the oil window, and using published oil viscosity-relax- 
ation charts to estimate viscosity. Another example is to 
use the map-derived bound fluid volume to compute 
Timur-Coates permeability according to the equation 
k_Timur = a * phit a b * ((phit bfv)/ bfv) Ac, where a, b, c 
are constants, phit is total porosity, and bfv is the bound 



fluid volume. 

[0060] a second approach to directly determine quantitative re- 
sults from D-T2 maps involves a determination of log- 
mean diffusion (D ). In standard MRF analysis, the raw 

LM 

data is fit directly using the constraints of Eqs. 1, 2 and 3. 
Similar constraints are imposed for other prior art model 
dependent inversions. An alternative approach, as dis- 
closed herein is to use the maps themselves as input to 
derive the solution which MRF attempts. Since the infor- 
mation contained in the maps is essentially identical to 
that of the original data, the two methods of solution 
should be comparable. In practice, however, the data is 
often lacking in diffusion information and therefore, the 
different fluids D-T2 amplitudes are spread over large ar- 
eas of the map (resolution). This is in contrast to the ideal 
situation for application of the point-and-click method 
where each fluid artifact is substantially separated from 
other fluid instances. The problem then consists of re- 
assigning the amplitude spread in the diffusion axis to the 
different formation fluids. According to this second ap- 
proach, an approximate way to do this is to use the geo- 
metric mean diffusion rate for each T2, notated D (T2), 

LM 

computed from the maps, and redistribute the amplitude 



at this T2 according to the chosen fluid model. For exam- 
ple, for a model consisting of water and oil it is conve- 
nient to define an apparent water saturation at each T2 
value, SXO(T2), 

£n« (T2)=D w (T2f* i " i DvlYl) 1 -** 0 *") (5) 

Separate water and oil T2 distributions, F , and F , 

K H20 OIL' 

can now be derived, 

Fn io (T2, ) = 3X0(12, p<F(T2, ) (7) 

( T2 . ) = < 1 -XX0(T2 l ))xF(T2, ) ^ 

F(T2;)=ZZ F ( T2 ^r T K.^ t ) (9) 

[0061] The D approach is most effective when the fluid model 

LM 

has only 2 components. In such a case, a best guess of 
the fluid model is provided to the software algorithm. In 
the case of a model involving more than 2 components, 
the extra components are successively eliminated from 
the D-T2 map. For example, for a water-oil-gas model, 
the gas component can be eliminated from the map using 
the Visual-Point-and-Click approach as described above 
and the D-T2 map re-normalized for the water-oil model. 
In practice, because the formation always contains irre- 
ducible (bound) water, the final reduced model is either 



water-oil or water-gas. 
[0062] FIG. 16 illustrates a graphical representation of an exem- 
plary D approach. Here, the D-T2 map shows two d if— 

LM 

ferent fluid instances. The fluid instance or artifact 304 is 
likely oil due to its proximity on the theoretical oil re- 
sponse overlay 306. According to one embodiment, quan- 
titative evaluation of fluid instance 304 can be determined 
directly from the D-T2 using the above described point- 
and-click method. In this case, after artifact 304 has been 
evaluated, it can be deleted and the map renormalized 
containing only the fluid instance 302. However, since ar- 
tifact 304 is not well separated from artifact 302, the re- 
sults of the point-and-click depend greatly on the delin- 
eation process. In another embodiment, the D lm approach 
may include both fluid instances 302 and 304. Turning 
attention to fluid instance 302, the analysis is not as clear. 
Specifically, the distribution associated with fluid instance 
302 lies between what one would expect for either water 
or oil, indicated by the water overlay 308 and the oil over- 
lay 306. Further, because the fluid type cannot be deter- 
mined, any integration over the distribution would not 
provide accurate formation evaluation answers. A solu- 
tion, according to one embodiment, is to compute the 



mean diffusion D 312 across the distribution 302. Each 

LM 

value of D is then used to re-assign the signal ampli- 

LM 

tude according to the separate fluid diffusion rates as in- 
dicated by Eq. 6. The re-assignment is linear in log space 
and based on the proximity with respect to the fluids dif- 
fusions responses. For example, if log(D ) = a, log(D )= 

wat oil 

b, and DLM = c, then the water saturation is S = (c b)/(a 

w 

b). 

[0063] FIG. 17 provides another example using the D approach 

LM 

derived for a D-T2 map distribution. Here again, the prior 
art T2 graph cannot resolve the fluid type of instance 302. 
Specifically, according to prior methods, the T2 graph at 
the bottom-left panel shows only the total curve 322 
(curves 324 and 326 are computed using the methods 
disclosed herein). Furthermore, the diffusion distribution 
graph in the upper-right panel similarly cannot resolve 
more than one fluid type as indicated by diffusion curve 
320. According to the disclosed D approach, the D is 

LM LM 

multiplied by the overall T2 distribution as per Eqs 7 and 
8 to yield the water distribution 324 and oil distribution 
326 seen in the bottom left panel. Integrations of the wa- 
ter and oil distributions then give the water volume (PhiW) 
and the oil volume (PhiO) seen in the bottom right pane. 



The oil log mean T2 can also be computed from the oil 
distribution (vertical dash line) from which the oil viscosity 
is estimated (Vis). 
[0064] FIG. 18 shows a flow chart of an exemplary method em- 
ploying a D determination. Beginning at step 800, a D- 

LM 

T2 map is generated or imported from a model indepen- 
dent inversion process. At step 802, a fluid model is se- 
lected by the user based in part on a visual inspection of 
the D-T2 map or prior knowledge of the sample, or a 
combination of both. The most common fluid models in- 
clude water-oil or water-gas, although other combina- 
tions are equally applicable. Additionally at step 802, any 
fluid artifacts not included in the selected fluid model, but 
present on the D-T2 map, are removed from the distribu- 
tion. This may either be performed manually using a com- 
puter mouse, for example, or automatically by software 
based on the proximity of an unwanted artifact to the the- 
oretical response of the fluid that is not part of the model. 
For example, for a selected water-oil model, software can 
detect that a concentration of amplitude, indicating the 
likely presence of a fluid, occurs near the theoretical gas 
response. An edge-detection routine is executed to deter- 
mine the extent of the artifact, interpreted as a gas, and 



then delete that region of the map from the total distribu- 
tion. 

[0065] Once the unwanted artifacts have been removed, the D-T2 
map at step 804 is normalized to the selected two-fluid 
model. At step 806, the log mean diffusion is calculated 
over the extent of the amplitude concentration. According 
to one embodiment, the D curve is displayed as an 

LM 

overlay on the D-T2 map. From there, and as described 
above, the fluid volume and saturation is determined at 
step 808 using Eq. 6 and 7. Continuing at step 810, the 
volume and fluid saturation indications are adjusted by 
focus windowing to improve the accuracy of the evalua- 
tion answers. Specifically, window focusing is performed 
to impose a saturation value, over a T2 region, in effect 
overriding the saturation computed by DLM. This is pre- 
ferred to counteract the unwanted effects of restricted 
diffusion, internal gradients etc. mentioned before. Fi- 
nally, at step 812, once fluid volume and saturation quan- 
tities, or other base values, have been determined, auxil- 
iary computations at 510 can be requested by the user or 
automatically by the algorithm. For example, at step 512, 
viscosity can be determined by computing the mean T2 of 
the oil distribution and using published oil viscosity-re- 



laxation charts to estimate viscosity. 

[0066] it should be noted that the above described exemplary 
approaches, including the point-and-click method and 
the D lm method may be employed on variations of the D- 
T2 maps. For example, it can be helpful to define a D-T2 
map for a certain T1/T2 ratio. In the context of the point- 
and-click method, the D-T2 inversion results are normally 
integrated over the third dimension which is T1/T2 ratio. 
However, it is also possible to decompose the total signal 
into separate results corresponding to each T1/T2 value 
and benefit from Tl information in the fluid interpreta- 
tion. Essentially, this means that a D-T2 map is generated 
for each T1/T2 plane. 

[0067] FIG. 19 provides an example where elimination of one 

fluid phase is necessary so that the D approach can be 

LM 

applied. The D-T2 map in the top left panel is derived 
from the D-T2 map shown in FIG. 14. The gas peak 240 
has been eliminated from the original map of Fig. 14 us- 
ing the visual point-and-click approach and the remaining 
amplitude is subsequently re-normalized to give the new 
map shown in FIG. 19. Specifically, the re-normalized map 
shows artifacts 242" and 244", both believed to be oils, 
and artifact 246", believed to be water but unclear. Note 



the absence of the gas diffusion amplitude 254 in the top 
right pane. Here again, neither diffusion peak 256" and 
258" fully resolve any of the artifacts. Specifically, each 
diffusion peak contains information from multiple arti- 
facts, peak 258" containing information from fluid in- 
stances 244" and 246", and to a lesser degree, peak 256" 
containing information from fluid instances 242" and 
246". Initial fluid volumes and saturations are displayed in 
the bottom right pane. However, from the D-T2 map and 
the T2 graph showing an estimated oil peak 248a and an 
estimated water peak 248b, it is apparent further inter- 
pretation is needed to better determine the fluid type as- 
sociated with artifact 246". 
[0068] | n FIG. 20, the peak appearing at short T2 (~ 5ms), peak 
248" in FIGS. 14 and 19 is believed to be water. The 
breadth of the peak in the D-T2 map reflects the uncer- 
tainty in the corresponding diffusion rate. This is due to 
lack of information in the original NMR data. With a priori 
knowledge we can manually assign this peak as water us- 
ing the window focusing method, resulting in peak 248"". 
Having done this, the water volume, oil volume and the T2 
log mean of oil have changed compared to their initial 
values shown in FIG. 19. 



[0069] FIG. 21 shows an example of interactive computation of 
the two oil viscosities that correspond the two bright 
spots along the oil response line 247 and 247" on D-T2 
map (FIGS. 14 and 19). Timur-Coates permeability is also 
computed interactively for a user-defined T2 cutoff. 

[0070] while the invention has been described with respect to a 
limited number of embodiments, those skilled in the art, 
having the benefit of this disclosure, will appreciate that 
other embodiments can be devised which do not depart 
from the scope of the invention as disclosed herein. For 
example, embodiments of the invention may be practiced 
with a wireline tool as well as a LWD or MWD tool. In addi- 
tion, embodiments of the invention may be practiced on a 
fluid sample removed by a formation tester and the NMR 
measurements are either acquired in the formation tester 
or in a laboratory. Further, the disclosed methods are not 
acquisition specific and may be applied to nearly all 
datasets regardless of whether a CPMG, diffusion editing 
or other pulse sequence is used. Accordingly, the scope of 
the invention should be limited only by the attached 
claims. 



